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ABSTRACT 



Wc present the first near infrared (NIR) spatially resolved images of the circumstcllar 
transitional disk around SR21. These images were obtained with the Subaru HiCIAO 



^Eureka Scientific, 2452 Delmer, Suite 100, Oakland CA 96002, USA 

^Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA 

'^H.L. Dodge Department of Physics and Astronomy, University of Oklahoma, 440 W Brooks St Norman, OK 
73019, USA 

^Department of Astrophysical Sciences, Princeton University, NJ 08544, USA 

®The Center for the Promotion of Integrated Sciences, The Graduate University for Advanced Studies (Sokendai), 
Shonan International Village, Hayama-cho, Miurargun, Kanagawa 240-0193, Japan 

^"Laboratoire Lagrange, UMR7293, Universite de Nice-Sophia Antipolis, CNRS, Observatoire de la Cote d'Azur, 

06300 Nice, France 

^'^Max Planck Institute for Astronomy, Konigstuhl 17, 69117 Heidelberg, Germany 

^■^Departinent of Physics and Astronomy, College of Charleston, 58 Coming St., Charleston, SC 29424, USA 

^^Department of Astronomy and Astrophysics, University of Toronto, 50 St. George Street M5S 3H4, Toronto 
Ontario Canada 

"Subaru Telescope, 650 North Aohoku Place, Hilo, HI 96720, USA 

^^Universitats-Sternwarte Munchen, Ludwig-Maximilians-Universitat, Scheinerstr. 1, 81679 Munchen, Germany 

^^Department of Astronomy, The University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan 

^''Institute for Astronomy, University of Hawaii, 640 North A'ohoku Place, Hilo, HI 96720, USA 

^^ExoPlanets and Stellar Astrophysics Laboratory, Code 667, Goddard Space Flight Center, Greenbelt, MD 20771 
USA 

^^Department of Astronomy, Kyoto University, Kitashirakawa^Oiwake-cho, Sakyo-ku, Kyoto, 606-8502, Japan 

^"OfRce of the President, Hiroshima University, 1-3-2 Kagamiyama, Higashi-Hiroshima, 739-8511, Japan 

^^Departamento de Astrofi'sica, CAB (INTA-CSIC), Instituto Nacional de Tecnica Aeroespacial, Torrejonde Ardoz, 
28850, Madrid, Spain 

^^Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, 91109 USA 

^^TMT Observatory Corporation, 1111 South Arroyo Parkway, Pasadena, CA 91105, USA 

^''institute of Astronomy and Astrophysics, Academia Sinica, P.O. Box 23-141, Taipei 106, Taiwan 

^^Astronomical Institute "Anton Pannekoek" , University of Amsterdam, Science Park 904, 1098 XH Amsterdam, 
The Netherlands 

^^Kavli Institute for the Physics and Mathematics of the Universe, The University of Tokyo, Kashiwa 277-8568, 
Japan 

^'^Department of Cosmosciences, Hokkaido University, Sapporo 060-0810, Japan 
^® Astronomical Institute, Tohoku University, Aoba, Sendai 980-8578, Japan 



- 3 - 



camera, adaptive optics and the polarized differential imaging (PDI) technique. We 
resolve the disk in scattered light at H-band for stellocentric C/l^r^O'.'G (12<r<75AU). 
We compare our results with previously published spatially-resolved 880//m continuum 
Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radia- 
tive transfer models reveal that the large disk depletion factor invoked to explain SR21's 
sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more mod- 
erate depletions (6=0.1, 0.01 relative to an undepleted disk) than those that reproduce 
the sub-mm cavity ((5~10~®) are inconsistent with our H-band images when they are as- 
sumed to carry over to small grains, suggesting that surface grains scattering in the NIR 
either survive or are generated by whatever mechanism is clearing the disk midplane. 
In fact, the radial polarized intensity profile of our H-band observations is smooth and 
steeply inwardly-increasing with no evidence of a break at the 36AU sub-mm 

cavity wall. We hypothesize that this profile is dominated by an optically thin disk 
envelope or atmosphere component. We also discuss the compatibility of our data with 
the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, 
and find that we can neither exclude nor verify this scenario. This study demonstrates 
the power of multiwavelength imaging of transitional disks to inform modeling efforts, 
including the debate over precisely what physical mechanism is responsible for clearing 
these disks of their large midplane grains. 



1. Introduction 



Transitional disks are of particular interest to the astronomical community because they have 
features intermediate between optically thick gas-rich protoplanetary disks and optically thin dusty 
debris disks, suggesting that they may be an intermediate evolutionary step. Transitional disks were 
first classified based on their Spectral Energy Distributions (SEDs) as objects with significant mid 
infrared excess emission indicative of starlight reprocessed by a circumstellar disk, but a lack of 
strong near infrared (NIR, ~l-5/im) excess. This NIR deficit implies a dearth of the very warmest 
grains in the inner disk (Strom et al. 1989: Skrutskie et al. 1990). Spatially-resolved sub-mm 



imaging of a number of transitional disks has revealed that many of them are indeed cleared of 
large thermally-emitting dust grains out to lOAU or more, well beyond the dust sublimation radius. 



thereby suggesting that an additional clearing mechanism may be at work (Andrews et al. 2009 



Hughes et al.||200"9t |Brown et al.||2009l [Andrews et al.|[2011[ referred to as All throughout the 
remainder of this paper). 

The precise mechanism by which these clearings are made is widely debated in the literature, 
with several viable alternatives put forward, including clearing by massive forming protoplanets 



^ Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of 
Japan. 
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(e.g. Dodson-Robinson &; Salyk 2011 Zhu et al. 2011), clearing by stellar mass companions (e.g., 



Ireland Sz Kraus 2008 Jensen & Mathieu 1997), photoevaporation (e.g. Owen et al. 2011, Clarke 



et al. 2001 Hollenbach et al. 1994) and grain growth (e.g., Birnstiel et al. 2012 Dullemond & 
Dominik 2005). Each theoretical explanation for the observed presence of sub- mm disk cavities 



makes its own predictions for clearing timescales, efficiencies and radii. However, the small amount 
of available spatially-resolved multiwavelength data on transitional disks has limited the ability of 
modelers to test these theories against observation. As such data accumulate, they reveal an in- 
creasingly complex picture of the "cleared" regions. In particular, cavities seen at long wavelengths 
do not necessarily manifest themselves in NIR scattered light (e.g. this work, Dong et al. 2012). 



Spectroscopic studies also reveal a variety of gaseous species in the inner regions of transitional 
disks (Salyk et al. 2009: Pontoppidan et ar]|2008 >. As large dust grains, small dust grains and gas 
all appear to respond quite differently to the mechanism or mechanisms that create transitional 
disks, a thorough understanding of their clearings can only be grasped through accumulation of a 
statistically significant sample of multiwavelength data on transitional disks. 

Thus far, the most success in spatially resolving transitional disks has been achieved in the 
sub- mm regime (e.g.. All), with a small number of additional studies using mid-infrared (MIR, 



10-40/im) aperture masking (e.g. [Eisner et al. 2009 hereafter E09), NIR spectroastrometry (e.g. 



Pontoppidan et al. 2008, hereafter P08) and MIR interferometry (e.g. Ratzka et al. 2007) to infer 



inner disk structure. Efforts are currently underway to image a sample of transitional disks at 
shorter wavelengths, including the Strategic Exploration of Exoplanets and Disks with Subaru 



(SEEDS; Tamura 2009) campaign. 



Despite initial success at resolving the rim of the LkCa 15 cavity in scattered light (Thai- 



mann et al. 2010) with a Locally Optimized Combination of Images (LOCI) method, additional 



transitional disks imaged by SEEDS are revealing that these cavities, including that of LkCa 15, 



are more elusive when viewed in polarized light (Hashimoto et al. 2012 Muto et al. 2012, Dong 



et al. 2012: Mayama et al. 2012 Grady et al. 2013, Wisniewski et al. 2013, in prep,). On the sur 



face, this suggests that NIR scattering material is present, and even abundant, inside the cavities. 
Although sub-mm studies are generally revealing that transitional objects are indeed cleared or 
depleted of large grain material inside of the cavity radii inferred from their SEDs, NIR 'clearings' 
in transitional disks may not be as ubiquitous nor as fully cleared as their sub-mm cousins. In fact, 
it appears that some transitional disks contain abundant small grain material inside of observed 
sub-mm cavities. 

Polarized NIR scattered light and sub-mm imaging, though both sensitive to circumstellar dust, 
probe very different regimes in the disk, both in terms of density and height above the midplane. 
The sub-mm emission from transitional disks is thermal, tracing mainly mm-sized grains that are 
concentrated ("settled") toward the disk midplane (e.g., Dullemond fc Dominik|[2004 Brauer et al 
2008t|D'Alessio et al.|200"6l ). These large (a>10/im) grains are very important for understanding the 
large scale structure of the disk, as they house most of the disk's dust mass and their coagulation 



drives the process of planet formation under the core accretion scenario (e.g.. Pollack et al. 1996) 
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The density and temperature of the grains at the disk midplane also determine the size and shape 
of the remainder of the disk. 

At NIR wavelengths, however, it is unclear whether scattered light emission originates from 
an optically thick small grain disk surface, or from more diffuse material in or above the disk. This 
uncertainty that is compounded by the small number of spatially resolved disks in the literature. It 
is unknown whether disks with large grain clearings in the sub-mm are also cleared of small grains, 
nor is it understood what proportion of the disk mass lies in the small grain population. The 
fact that we see axisymmetric polarization structure in the NIR from transitional objects reveals 



that they house circumstellar disks with some small grain material (Whitney & Hartmann 1992). 
However, very little dust mass is needed to scatter at H-band, so it is unclear how much mass resides 
in this portion of the disk and whether it is optically thin or optically thick in this wavelength regime. 
If transitional disks are optically thick in the NIR, as is the case for protoplanetary disks, then 
this scattered emission will probe only the r=l surface. Small grains are believed to be vertically 



extended by a factor of two or more relative to larger grains ( Dullemond &: Dominik||2004 Brauer 
et al.||2008D , so NIR scattering in an optically thick disk should occur well above the disk midplane. 
It is unclear whether or not clearing mechanisms acting on the large grains at the disk midplane 
would have an observable effect on the NIR emission under this scenario. 

Conversely, if transitional disks are optically thin in this wavelength regime, as is the case for 
debris disks, then scattered light emission probes the entire vertical extent of the disk. In this 
case, scattered light emission originates from small grains throughout the vertical extent of the 
disk, down to and including the disk midplane. Observations of the same disk at both sub-mm and 
NIR wavelengths are compelling in that they allow for comparison of the settled and upper surface 
layers of the disk, the large (mm-sized) and small (/xm-sized) grains, separately. A good disk model 
should thus be able to self-consistently reproduce observations at both wavelengths simultaneously, 
as well as any additional multiwavelength data. Our efforts to do just that for the SR21 disk are 
described in detail in Section 4. Section 2 describes our data collection methodology, and our result 
are summarized in Section 3. 



1.1. Observational History of SR21 



SR21A is a 2.5M0, ~1 Myr old Class II Young Stellar Object (YSO) in the Ophiucus Star 
Forming Region. Recent astrometric observations suggest that the distance to the Ophiuchus 
star forming region is 119ib6 pc (Lombardi et al. 2008). However, we adopt a distance of 135pc 



( Mamajek||2008[ ) in order to remain consistent with previous studies. SR21 has been classified as a 



6'.'7 binary (Barsony et al. 2003); however, Prato et al. (2003) concluded that the two components 



were not coeval, and that the B component does not exhibit infrared excess emission and is therefore 
unlikely to host a circumstellar disk of its own. The A and B components also have discrepant 
proper motions ( Roeser et al.|[2010 ), so the "binary" does not appear to be a bound, coeval system. 
We will therefore refer to SR21A as simply SR21 throughout the remainder of this paper. 
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SR21 was first identified as a transitional disk by Brown et al. (2007), who inferred a 0.45-18AU 
gap in the disk based on a fit to the ~l-100/um SED. In order to explain the NIR excess seen at 
<5//m, they invoked a small uncleared region in the inner disk from 0.25-0.45AU. Since this initial 
identification, a variety of puzzling and seemingly contradictory data have emerged that continue 
to favor a disk that is not simply cleared of all material in this 0.45-18AU region. 

The first evidence of a lack of complete clearing in the inner disk came from CO rovibrational 
spectroastrometry done with CRIRES on the VLT by P08, who observed molecular gas emission 
originating from a region well inside of the dust gap inferred by Brown et al. They concluded that 
this emission arose from a narrow 7-7. 5AU ring of molecular gas. It is not surprising to have found 
gas emission inside of the dust gap, as SR21 may still be accreting, albeit weakly (<lO~^'^^M0/yr, 
a relatively high upper limit, Natta et al. 2006). However, the physical mechanism causing an 



apparently truncated ring of molecular gas at this radius is still an open question. P08 suggest that 
the inner 7AU cutoff corresponds to a physical truncation of the disk, while a rapid drop in gas 
temperature outward of 7AU, perhaps due to a thick disk component in the region, causes the CO 
flux to quickly drop below the detection threshold. This explains the observed truncation of CO 
emission at 7.5AU without requiring that the disk material be physically truncated at this radius. 

E09 also postulated that there is emission from within the ~18AU 'cavity' based on MIR (8.8 
and ll.G^m) aperture masking data taken with TReCS on Gemini South. They put forward an 



alternate interpretation to Brown et al. (2007), capable of fitting both the SR21 SED and their 



visibility data. They suggest that the disk that is completely cleared within lOAU, but with a 
warm companion at 10-18AU whose circum(sub)stellar material creates the NIR excess seen in the 
SED. 

The first spatially resolved imaging study of SR21 was done at 340GHz (880/im) and was 



published by Brown et al. ( 2009 ) among a sample of four disks with large central clearings. Followup 



data, also at 340GHz, were presented in Andrews et al. (2009). We base our sub- mm modeling 



efforts on All, a detailed imaging and modeling study of 12 transitional disks, including SR21. 
Using two-dimensional Monte Carlo radiative transfer modeling to fit the sub-mm data and the 
infrared SED, All concluded that the 36AU sub-mm central cavity in SR21 was depleted in material 
by a factor of ~10~® relative to the outer disk. 

As multiwavelength data accumulate on SR21, it is clear that its structure, particularly in the 
region interior to the 36AU sub-mm clearing, cannot be understood through SED modeling alone 



because there are several degenerate solutions capable of explaining it (Brown et al. 2007, Eisner 



et al.|2009 ). Spatially resolved images are needed at multiple wavelengths, probing multiple regions 
of the disk, before good theoretical constraints can be put on the mechanism clearing the interior 
of SR21. In the next section, we describe new spatially resolved NIR scattered light imaging, which 
contributes further to the puzzle of the SR21 "cavity". 
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Observations and Data Reduction 



Polarized Differential Imaging (PDI) of SR21 was done in H-band (1.6/xm) on 2011 May 22 
with the high-contrast imager HiCIAO (Tamura et al. 2006; Hodapp et al. 2006| ) and adaptive 



optics system A0188 (Minowa et al. 2010) on the 8.2m Subaru telescope on Mauna Kea. These 
observations were conducted as part of the larger SEEDS survey, which was begun in October 2009 
and is currently entering its third year. 

The data were taken in "qPDI" mode in which a double Wollaston prism is used to split 
the beam into four 512x512 pixel channels (O'.'0095/pixel), two each of o and e-polarizations. The 
splitting of each polarization state into two separate channels reduces saturation effects on the 
Hawaii 2RG detector. In order to obtain full polarization coverage and minimize artifacts, a half 
waveplate was rotated to four different angular positions (0°, 45°, 22?5 and 67?5). This cycle was 
repeated 18 times on SR21, with a 15 second exposure per waveplate position, for a total of 18 
minutes of integration time. 

Each image was bias subtracted, flat fielded and bad pixels removed in the standard manner for 
infrared data analysis using custom IDL scripts. Pinhole data were then used to create a distortion 
solution for each channel, which was applied before the images were rotated to a common on-sky 
geometry. All of the 15 second science exposures of SR21 were saturated to a radius of 6-8 pixels. 
For this reason, all subsequent combinations have had an inner r=8 pixel mask applied. 

Each individual channel was then aligned using Fourier cross-correlation. These aligned images 
were combined to create Stokes Q and U images using standard differential polarimetry methods, 
namely adding together each set of two identically polarized channels and then subtracting these 
combinations from one another. The four channels were also added directly together to form a 
total intensity image. 

To estimate the performance of the Subaru AO system during our observations, we used two 
unsaturated images of the Point Spread Function (PSF) reference star HD 148212, and 12 and 
10 unsaturated images of SR21 taken before and after the main SR21 data set respectively. The 
average Strehl Ratio of the HD 148212 images is 0.27, and the average Strehl Ratios of the SR21 
before and after images are 0.22 and 0.20 respectively. The standard deviations of the Strehl 
Ratios in the two SR21 data sets are 0.018 and 0.016 respectively, indicating that AO correction 
was maintained steadily. 



2.1. Polarized Intensity Images and Polarization Vectors 

The polarized intensity image shown in Figure 1 is a median combination of individual Po- 
larized Intensity (PI) images for each of the 18 waveplate cycles, computed via 
Polarization angles were calculated for each image according to the formula 9p = 0.5tan''^{U/Q) 
and median combined. These are shown overplotted on the PI image in five pixel bins in the 
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lefthand panel of Figure 2. The length of each vector has been scaled according to the magnitude 
of the polarized intensity in that bin. A first order instrumental polarization correction has been 



applied, as described in detail in Hashimoto et al. (2011). 



As has been demonstrated by other groups ( Close et al. 1997 Oppenheimer et al. 2008 Apai 



et al. 2004 Quanz et al. 2011), PI images are a powerful data product for studying circumstellar 
disks. Unlike most other varieties of disk observations, PI data are generally spared the difficult and 
imperfect task of PSF subtraction under the assumption that direct starlight is randomly polarized, 
while light scattered off of a circumstellar disk shows a preferred polarization direction according 
to the scattering location in the disk. Specifically, the scattering process will preferentially scatter 
light polarized perpendicular to the direction to the illuminating source, while light coming directly 
from the star will be unpolarized and randomly oriented. Extended polarized emission is therefore 
unlikely to occur in systems without scattered light disks, and light scattered from a disk should 
show a centrosymmetric pattern of polarization vectors. 

Although direct starlight is randomly polarized when it reaches the telescope, multiple reflective 
surfaces in the HiCIAO optical path are capable of polarizing some fraction of this incoming light. 
This phenomenon is referred to as instrumental polarization (IP). In order to estimate the quality 
of IP removal via standard techniques in our data, we reduced our data on the PSF star HD 148212 
using the same methodology as SR21, including a first-order instrumental polarization correction. 
HD 148212 is a G2V star at d=105±15pc and ~10' from SR21. It was observed just before SR21 
on the same night and in the same observing mode, with 17 cycles of waveplate rotations and 5 
seconds per exposure, for a total of 5.7 minutes of exposure time. As shown in the righthand panel 
of Figure 2, no polarized emission was observed beyond r>Cf!2 in the initial first order IP corrected 
HD 148212 polarized intensity image. Although the polarization signal inside of this radius is small 
compared to the Polarized Intensity of SR21 in the same region, it does show the characteristic 
four-lobed astigmatic pattern of the HiCIAO PSF. Additionally, while polarization vectors for HD 
148212 do not show centrosymmetry, they do show some alignment along the vertical axis of the 
instrument. These are both hints that there is a residual uncorrected instrumental polarization 
signature, and further correction is warranted. 

To first order, we can approximate the residual instrumental polarization by measuring the 
average polarization strength and direction of the signal in the HD 148212 polarized intensity 
image. We assume that HD 148212 is an ordinary diskless unpolarized star, and that all of its 
extended emission is due to residual instrumental polarization or photon noise. We construct 
an artificial "halo" by multiplying the average strength (1.1%) in the HD 148212 halo into the 
raw intensity image, then creating artificial "halo Q" and "halo U" images from it using the 
measured average polarization angle (-11°). We combine these synthetic Q and U images to create 
a "halo PI" image and subtract it from the original HD 148212 PI image. We also subtract the 
"halo Q" and "halo U" images from the raw Q and U images, and use the resulting images to 
calculate corrected polarization vectors, which we overplot on the halo-subr acted PI image for HD 
148212 in the righthand panel of Figure 2 . The extended polarized intensity emission from the 
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PSF star is reduced to r<G/!l through this process, and the resulting vectors are randomized, no 
longer showing a preferred polarization direction. We are therefore confident that this provides a 
reasonable approximation of the residual instrumental polarization. 

The need to apply this same "halo removal" procedure to the SR21 PI data can be seen in 
the Figure 2 SR21 polarization map. Although the SR21 vector pattern is clearly ccntrosymmetric 
along and around the disk major axis even in the uncorrected image, it deviates somewhat from 
centrosymmetry along the same instrumental axis as the aligned vectors in the HD 148212 image 
(approximately the minor axis of the disk). Since not all SR21 polarized emission is spurious, 
approximating the strength and direction of the halo is nontrivial in this case. We used a region 
with low intrinsic disk emission to construct a "halo" for subtraction, namely 100 pixel (0'.'95) wide 
regions centered on the minor axis at 30 pixels<r<50 pixels (0'.'285<r<0'.'475) on either side of the 
disk. Following the same procedure as with HD 148212, we created polarized intensity "halo" Q, U 
and PI images from the polarization direction (-8.8°) and strength (1.02%) measured in this region 
and subtracted. 

Subtraction of this first-order approximation reveals a more centrosymmetric pattern, as ev- 
idenced in the lefthand panel of Figure 3. However, it also introduces significant PSF artifacts 
to the polarized intensity map, as can be seen in the background polarized intensity contours of 
Figure 3. For this reason, we have chosen to model the "raw" unsubtracted radial profile in Section 
4, noting, however, that some of the emission interior to O'.'IS is likely to be contaminated by a 
small degree of instrumentally polarized stellar flux. The polarized halo strength of 1.02% is an 
order of magnitude smaller than our lower limit for the intrinsic SR21 polarization (10%), so this 
contamination is minimal. 

Despite its crudeness, the "halo subtraction" method highlights an interesting feature of the 
SR21 polarization vector structure, which is that several regions in the disk with high signal (notably 
the northeast and southwest quadrants) show deviations from centrosymmetry even after halo 
subtraction. This may indicate that there is unresolved structure in the disk at these locations. 
Interestingly, the deviation in the northeast is consistent with the location of the companion inferred 
by E09, which could be explained if the companion had an accretion disk of its own, with its own 
polarization structure superimposed on that of the SR21 disk. 

2.2. Intensity Images 

Intensity images were also constructed from the PDI observations of SR21 by cross correlating, 
centroiding and aligning the images in each of the four channels at each of the four waveplate (0° , 
45°, 22?5 and 67?5) positions for each of the 18 observation cycles. The four channels at each 
waveplate position were summed for each cycle, and the resulting 18 images were then median 
combined to create an SR21 intensity image. 

Unlike polarized intensity images, total intensity images are dominated by stellar flux, and 
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PSF subtraction is necessary. For this purpose, we created intensity images of HD 148212 using 
the same methodology described above. The final reduced (but unsubtracted) median combinations 
of the SR21 and HD 148212 total intensity images are shown in Figure 4. Both reveal significant 
uncorrected astigmatism, as evidenced by the four-lobed pattern centered on the star. The SR21 
disk is not visible in the unsubtracted total intensity image. The PSF subtraction process is 
described in detail in Section 3.3. 



3. Results and Analysis 

Our principle observational results are the final polarized intensity (Figure 1) and total intensity 
(Figure 4) images and polarization vectors (Figures 2 and 3). These images were analyzed in order 
to estimate the extent, inclination, major axis orientation and polarization percentage of the SR21 
disk. The derivation of each of these properties is described in detail in this section. 

3.1. Polarized Intensity Isophotes 

We began by fitting elliptical isophotes to our PI data in order to provide our own measure of 
the inclination and major axis of the SR21 disk, independent of the values in the literature, which 
were determined at other wavelengths. These elliptical fits are shown overplotted on Figure 5 in 
white, beginning at a radius of 10 pixels from the center (well outside the saturation radius) and 
extending outward to 50 pixels, where the PI flux approaches the noise level. Errors in isophotal 
fits were determined by summing in quadrature the intrinsic error estimation in the IRAFj^ ellipse 
procedure, which was used to fit the isophotes, and the standard deviation of the parameter across 
the image. 

The average ellipticity of these isophotes suggests a disk inclination of 14zb2° (where 0° would 
be a face-on geometry). This is smaller than the value reported by All; however, their estimate 
of 22° was made with the caveat that the sub-mm ring seen in SR21 is very narrow, and the 
geometry is not well constrained by their observations. The H-band emission is significantly more 
spatially extended; therefore, we adopt 14° as the preferred value in our modeling efforts, which are 
described in section 4. There is no statistically significant deviation in our calculated inclination 
value between the inner and outer disk of SR21. 

There is a clear counterclockwise progression from the inner to the outer disk in the major 
axis of the elliptical isophotal fits (overplotted in white in Figure 5). The average orientation of 
the major axis over the entire disk is 86ibll°. The large uncertainty is due to the fact that the PA 
migrates from inner to outer disk. This happens relatively smoothly, but with a pronounced jump 



^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy (AURA) under cooperative agreement with the National Science Foundation. 
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between 17 and 20 pixels. Inside of this radius, which is of particular interest because it corresponds 
to the location of E09's hypothesized companion, the average major axis position angle is 76ib8° 
(measured East of North), while outside of this radius it is 96ib7°. This shift in the average major 
axis PA of ~20° between these two regions, and of 42° in total between inner (68°) and outermost 
(110°) isophotes, suggests that the SR21 disk may contain unresolved structures such as a warp or 
a spiral arm. Some of this PA migration might be explained by instrument ally polarized emission 
with the shape of the HiCIAO PSF, particularly because the most pronounced turnover occurs at 
approximately the same radius as the extent of the HD 148212 halo. However, the observed SR21 
isophote major axis migration continues well beyond that radius, and we therefore conclude that 
it is a real feature of the data. 

We combine our results with previous estimates of the major axis PA from the literature in 
Table 1. The All estimate for the major axis PA of 100° (overplotted on Figure 5 in red) is 
consistent with our measurement of the PA of the outer disk. We find that previous estimates of 
the PA of the disk are strongly dependent on the radius being probed, and our observation of a 
counterclockwise trend with radius fits well into the broader picture. This suggests that there is 
unresolved non-axisymmetric structure in the SR21 disk, and provides further confirmation that 
the trend we observe is not caused purely by residual PSF-shaped instrumental polarization signal. 



3.2. Radial Brightness Profiles and Asymmetries 

We used the final median combination of unsaturated total intensity images and the published 
2MASS photometry values of SR21 and HD 148212 to establish a flux conversion for our data. In 
doing so, we make the assumption that the H-band magnitudes of both objects remained steady in 
the decade between 2MASS and HiCIAO observations. H-band emission of T Tauri stars has been 



shown to be variable (Carpenter et al. 2002), however the The ADU-to-Jansky conversion factor 
calculated using the HD 148212 images differs by only 1.6% from that calculated using the SR21 
images; therefore, we believe that this is a robust assumption and have adopted the average of the 
two values as our flux conversion factor. We have applied a software mask to the saturated region 
interior to r=8 pixels; however, astigmatic PSF artifacts extend out to r~13 pixels (O'.'IS, ~18AU). 
Our data are therefore most robust at 18AU<r<80AU. 

Computed flux-converted radial polarized intensity profiles for SR21 are shown in Figure 6. 
Brightness profiles are shown for both the azimuthally-averaged case (filled black points) and for 
the case where we limit analysis to a region 11 pixels wide centered on the disk major axis (open 
red points). We also over plot (open blue points) the radial profile of a PSF-subtracted Polarized 
Intensity image. Although PSF subtraction is nominally unnecessary for polarized intensity images, 
our analysis in section 2.1 revealed that residual instrumental polarization contributes at a low level. 
We therefore utilized the polarized intensity image of the PSF star HD 148212, and scaled it such 
that the noise in the background region was equivalent to that of SR21 and subtracted. This 
provided a subtraction with fewer evident PSF artifacts than our earlier "halo subtracted" PI 
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image and serves to verify that the slope and steepness of the radial polarized intensity profile is 
insensitive to contamination by residual instrumental polarization. 

In all three cases, we plotted the profile in five pixel radial bins excluding the innermost 

two bins (so interior to CZ/l), which have pixels inside the saturation radius. This binning was 
chosen to match the diffraction limit of (^!05 at H-band. The brightness profile along the major 
axis is systematically higher than the azimuthally-averaged profile, as is to be expected based on 
the increased brightness of the disk in this region; however, the shapes are virtually identical. 
The PSF subtracted profile is systematically lower than that of the unsubtracted profile, as is 
to be expected. Unlike the other two profiles, the PSF-subtracted profile does not flatten in the 
outer disk, suggesting that this trend may be caused by residual instrumental polarization in the 
unsubtracted profiles. We have elected to fit the azimuthally-averaged profile of the unsubtracted 
image in our modeling, which is described in Section 4. 

These radial PI profiles reveal that H-band scattered light emission exists in the SR21 disk 
out to a distance of at least ~(y.'6, which is comparable to the disk extent in the All SMA images. 
This corresponds to a physical scale of ~80AU, or more than twice the radius of the observed 
sub-mm cavity. The excess in the NE quadrant of the disk, which could be due to the contribution 
of a companion as discussed in Section 3.3, is outside of the region used to create the major axis 
radial profile, and is smoothed out in the azimuthally averaged profile, therefore does not have any 
noticeable effect on either of our derived radial profiles. 

Errors in the radial polarized intensity measurements were estimated using both the back- 
ground photon noise in the four individual qPDI channels of the SR21 raw images and using the 
observations of the PSF star HD 148212. The residual polarized intensity flux surrounding HD 
148212 was summed in the same manner as the SR21 radial profile and is comparable to the pho- 
ton noise. The error bars in Figure 6 and all of the radial profile plots shown in the modeling 
section correspond to the photon noise. 

A canonical power law for the surface brightness of a flared disk is shown overplotted in 
green on all radial proflles and is a poor flt to the data. The best flt radial profile to the data is 
of the form r~^. Although there is some indication that it may fiatten beyond 0'.'45, this trend is 
not exhibited in the PSF subtracted PI profile, so we believe it to be spurious. The smooth, steep 
incline in polarized intensity from the outer disk in to the saturation radius is the crux of the new 
NIR data, and explaining it is the main driver of the modeling efforts described in Section 4. 

3.3. Intensity Images and Compatibility with the Existence of a Substellar 

Companion 

Direct addition of the 4 qPDI channels yields a total intensity image, which may also be 
diagnostic of the disk, but is dominated by stellar emission. PSF subtraction is required in order to 
isolate the disk component in these images. PSF subtraction is an imperfect process in systems such 
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as SR21 for a number of reasons. First, Strehl Ratio may vary between the two sets of observations, 
and within each set itself, and Strehl Ratio mismatch can introduce artificial structure into the PSF- 
subtraction. The SR21 and HD 148212 data sets were calculated to have comparable Strehl Ratios 
(~20%); therefore, we do not believe that Strehl mismatch is a dominant effect. Secondly, speckles 
that survive the median combination process are not necessarily cospatial in the two data sets; 
therefore, local excesses and deficits must be viewed with some suspicion. Thirdly, the SR21 data 
set was centrally saturated, while the HD 148212 data set was not, and so the counts in the SR21 
image are an underestimate of the total emission. Finally, and perhaps most importantly, the total 
intensity of SR21 includes both a direct stellar and scattered light disk contribution, while tHD 
148212's emission corresponds to direct starlight alone. These last two effects make estimation 
of the appropriate scaling of the PSF image for subtraction somewhat ad hoc, but by employing 
a series of reasonable approximations we can place bounds on the appropriate scale factor for 
subtraction. 

We chose to scale the PSF star HD 148212 first by the relative ratio between its 2MASS 
H-band magnitude and that of SR21, a quantity which we will call a. Assuming that the total H- 
band magnitude of these systems has not varied significantly with time, subtracting the PSF star 
intensity image scaled by this value is an oversubtraction, both because some portion of SR21's 
H-band magnitude corresponds to its disk contribution and because the SR21 image is saturated 
at the center. Therefore we use this subtraction (SR21 - a*HD 148212) as an "upper limit" on the 
disk's total intensity in scattered light. 

We then created a grid of PSF subtractions with successively lower scale factors applied to 
the PSF star. We used the ratio of scattered polarized light to total scattered light in the disk as 
a metric in establishing the lower limit for the PSF scale factor. This ratio of polarized to total 
scattered light (P=Pl/ldisk) should be on the order of 50% in the case of 100% efficient scattering 
of unpolarized starlight off of spherical grains; however, variables such as porosity, oblateness, and 
alignment as well as particle size and composition effects, multiple scattering effects, wavelength 
dependence and inclination effects have been invoked to explain observed deviations from this value 



(Whitney & Wolff 2002, Simpson et al. 2009, e.g). Observed scattered polarization fractions vary 
from ~10% to ~70% in transitional disks, with typical values of approximately 40% (Dean Hines 
and Glenn Schneider, private communication based on HST program GO 10178, October 2012). 

We deem values of <10% to be unphysically low based on these results; therefore, we apply 
the scale factor resulting in P=10% (0.5a) as our lower limit for PSF subtraction. Shown in Figure 
7 is the PSF subtraction resulting from a scale factor 0.75a, which is halfway between the upper 
(a) and lower (0.5a) limits established above. The image has been smoothed by a gaussian with 
FWHM=6 pixels to mitigate the effects of speckles on the inferred structure of the disk. 

In estimating the ratio of polarized to total disk intensity P (PI /Idisk), we use this intermediate 
PSF scaling and an annulus of 15<r<60 pixels. This region begins just outside the radius of the 
PSF astigmatism and ends at approximately the outer radius of the polarized intensity disk. Within 
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this region, we have rejected all pixels where (i) the polarized intensity flux is below noise level, (ii) 
P>1.0, indicating a PI speckle, or (iii) P<0, indicating an I speckle. Our upper and lower limits 
for the subtraction tell us that SR21's P value lies somewhere between 10 and 27%, and a=0.75 
gives P=18%. 

The final PSF-subtracted and masked intensity images of SR21 for all scale factors investigated 
reveal two structures, both of which are also present in the PI data. The first is an excess in the 
NE quadrant near the location of the E09 companion, and the second is an arc extending upwards 
towards the north. These are marked as A and B in Figure 7. Although it is entirely possible that 
these features are spurious, they are intriguing in the context of the E09 observations and our own 
observation of the trend in PA with radius in the SR21 disk, which is suggestive of structure. 

E09 inferred a warm (~700 K), extended (~40 Rq) companion 18 AU from SR21 based on the 
system's mid-IR flux and visibilities. Such a source would contribute ~5% of SR21's bolometric 
flux, or ~0.4 Lq, and would need to be a low- mass star surrounded by a dusty shroud or disk. 
The enshrouded companion would likely be considerably brighter in the i?-band than a ~700 K 
blackbody, itself a magnitude 16 source at SR21's distance. We note that a bolometric luminosity 
of 0.4 Lq is characteristic of a ~0.2 Mq T-Tauri star at 1 Myr, which would have an effective 
temperature of ~3000 K ( Kenyon Sz IIartmann|[l995 ) . The observed excess at point A in our data 



could also be due to a simple overdensity in this region of the disk. 

We obtained our images of SR21 with the image rotator off, allowing us to use Angular Differ- 



ential Imaging (ADI) to search for a companion. A traditional ADI analysis (Marois et al. 2006) 
detected no statistically significant sources. Because of SR21's declination, we were only able to 
obtain ~10° of field rotation, making it difficult to apply more sensitive techniques like Locally 



Optimized Combination of Images (LOCI, [Lafreniere et al.|[2007 ). We therefore attempted a varia- 



tion of LOCI using the ACORNS- ADI pipeline (Brandt et al. 2012), in which we supplmented the 



comparison frames with 30 principal component images culled from a large library of observations 
of single stars. 

We detect no statistically significant companion candidates down to a limiting contrast of 100 
at 0//14 (18 AU at the distance to SR21). This is not a meaningful limit for a substellar companion. 
However, it does likely rule out a stellar companion sufficiently luminous to account for the mid-IR 
excess reported by E09. Follow-up observations at longer wavelengths or with extreme adaptive 
optics will be necessary to reach fainter detection limits. 



4. Modeling 



To assist in interpreting our results, we ran a suite of disk models using the three dimensional 
radiative transfer Whitney code (see Whitney et al. 2003a|b for details). The key capabilities, 
modifications from previously reported versions and input assumptions of the Whitney code are 
described in Section 4.1 
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In Section 4.2, we describe our efforts to reproduce the sub-mm observations of All using an 
independent code, and our verification of the All best fit parameters for the large grain dust disk. 

In Section 4.3, we report on several degenerate explanations for the three features that make 
the SED of SR21 unique: (1) the steep 12-20^m rise, (2) the marked NIR excess and (3) the very 
small size of the 10/im silicate feature. 

In Section 4.4 we report on our attempts to reproduce the observed steep polarized intensity 
radial profile, and its ability to constrain the geometry and optical depth of the disk. 

Section 4.5 details our efforts to merge these three pieces into a single self-consistent model, 
and in Section 5 (Discussion) we speculate on the implications of our best fit model. 



4.1. The Whitney Code and Input Assumptions 



The Whitney code mimics settling and grain growth by specifying the composition and struc- 
ture of a "large grain" disk separately from that of a "small grain" disk. The large grain disk can 
be artificially settled relative to the small grain disk by setting its scale height parameters to lower 
values. Although both disks contain a distribution of grain sizes, we will refer to them simply as 
the "large grain disk" and "small grain disk" throughout. 

The Whitney code treats both thermal emission and scattered light in the disk. Grain tem- 
peratures are computed using the method of |Lucy||1999 and the code treats anisotropic scattering 



using a Henyey-Greenstein function (Henyey &; Greenstein 1941). The albedo, average scattering 



angle and maximum polarization are pre-computed at each wavelength based on a Mie calculation 
for the ensemble of grains in each population. The polarization phase function is approximately 
dipolar, following |White|[T979 



The Whitney code utilized herein has been modified from previously reported versions and is 



described in detail in Dong et al. (2012) and Whitney et al. 2013 (in prep). The modifications of 



particular relevance to this paper were threefold. First, the disk density distribution was modified 
in order to allow for an exponential cutoff in mass as a function of radius, making it identical to 
the disk mass distribution of models run by All. Secondly, the new models allow for a "puffing 
up" of the inner disk edge and the inner edge of the disk gap, a physical situation widely invoked to 
explain NIR spectral excess (e.g., Dullemond Sz MonnierpOlO ) and increasingly invoked to explain 



the observed features of gapped disks (e.g., Jang-Condell Sz Turner 2012). The final modification 
made to the code allows for the differential depletion of large and small grains, in addition to 
independent geometries between the two disks (a feature included in earlier versions of the code). 
This means that the disk radial density exponent (a, pocr~"), radial scale height exponent (/3, 
Z(xr~^), absolute scale height at lOOAU and depletion (5) are all specified independently inside 
and outside of a disk gap or cavity. In effect, this allows for four independent disk components - 
large and small grain disks inside and outside of a specified gap radius. The size distribution and 
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composition (and therefore opacity as a function of wavelength) of large and small grains are also 
specified independently. 

The basic model input parameters (e.g., stellar mass, radius, temperature) were determined 
based on values in the literature, and those adopted by All in particular. The choices of model 
parameters are given in Table 2. In total, over 250 independent models of the SR21 disk were 
calculated in this modeling effort. 

Initial input assumptions included a disk mass of O.OOGMq, a gas:dust ratio of 100:1, and a 
large grain dust mass fraction of 85%. Neither of these values is well constrained by our H-band 
observations, so we have assumed the same values as All. We have fixed the composition of the 
settled dust disk to the grain prescription of ( Wood et al.|2002 ), a power-law grain size distribution 



with an exponential cutoff (Kim et al. 1994) composed of amorphous carbon and astronomical 
silicates with solar abundance constraints for carbon and silicon (C/H ~ 3.5 xlO^ and Si/H ~ 3.6 
xlO^, respectively; Anders & Grevesse 1989 Grevesse & Noels 1993). The functional form of this 
size distribution is 



n{a)da = Cia ^exp 



da 



, where a is the grain size, n is the number of grains, and the parameters p, Oc, and q, which 
control the distribution shape, are adjusted to fit the wavelength dependence of the observations. 
In particular, the long wavelength slope of the SED and the peak of the MIR excess are sensitive to 
these parameters. Following the [Wood et al.|[2002 best fit model for the HH30 disk, we use p=3.0, 
q=0.6 and ac=50/xm. This is a good match to the grain properties inferred by All. 

Several small grain dust prescriptions were also investigated, and the observed SED and H- 
band images were found to be relatively insensitive to their particulars. We use the interstellar 
medium dust prescription of Kim et ar]|1994 which has a smaller average grain size than our other 
models. Like the large grain disk, the small grains follow a power law size distribution with an 
exponential cutoff, however in this case they are composed of silicates and graphite, with maximum 
grain sizes of 0.16/im and 0.28;um respectively. 



4.2. Sub-mm Image Modeling 

Our H-band PDI observations probe only scattered light from the upper surface of the disk 
and therefore do not constrain the properties of the large grain disk where the majority of the disk 
mass lies. However, we have endeavored to match the sub-mm observations of All in addition to 
our newly- acquired data in order to develop the most physically motivated disk model possible. As 
a first test of the compatibility of the All models with the Whitney code, we began with input 
parameters for the disk that were identical (or nearly) to the best-fit model of All. 

Figure 8 shows the SMA image taken by All in panel 1, and the equivalent 880//m image 
produced by our model in panel 2. We are able to produce a similarly well- fitting model with 
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the Whitney code and identical input parameters. Model outputs were sampled with the same 
spatial frequencies observed with the SMA, and images were created with the same inversion and 
deconvolution parameters as the observations. 

Having shown that the All and Whitney models produce equivalent results with identical 
input assumptions, we then investigated the effect of gap depletion on the sub- mm images to test 
the robustness of the All assumption that the gap is depleted by a factor of 6=10~^'^ relative to 
the outer disk. Keeping all other disk parameters fixed, we varied only the depletion factor 6 in the 
interior of the gap by an order of magnitude for each model from (5=10~^ (two orders of magnitude 
higher depletion than the All assumption of ~10~^) to zero depletion (which, at least on the 
surface, is the most consistent with the evidently undepleted cavity seen in our observations). 

We break these models into three categories: "high" depletion models (5=10^^-10~^), "mod- 
erate" depletion models ((5=10~^-10~^) and "low" depletion models ((5=10~^-10°) and show the 
real 880/im image of All, our simulated image, residuals and visibility fits for each of these models 
in Figures 9-11. In each case, the model images were scaled such that the total flux was equivalent 
to that of the real SMA observation, and the required scaling is given in the model panel of the 
figures. 

The "low" depletion models provide poor fits to both the 880/xni image and the visibility 
profile. We conclude, as did All, that the SR21 disk is depleted in large sub-mm emitting grains 
interior to 36AU. All of the moderate and high depletion models provide reasonable fits to both 
the 880//m image and the visibility profile; however, the S=10~^ model shows the lowest residuals 
and the 5=10^^ model provides the best fit to the visibility profile. The 5=10^^ model requires 
just a 2% increase in flux to sum to the observed flux density (easily explained by small differences 
in input opacity assumptions), while the 5=10~^ model requires an 8% scaling, so we adopt the 
S=10~^ value for the large grain r<36AU disk in all further modeling. This value is closest to the 
All best fit value of d=10~^-^; therefore, we consider our large grain models to be functionally 
equivalent. 

As a test of the robustness of the large grain disk assumptions, we verified that the small 
grain disk parameters specified in our models did not create observable variation in the simulated 
880/7,m image. As expected, we found that wide variation in input small grain parameters has no 
observable effect on the 880/im images produced by our models. The 880/xm emission of the large 
grain disk is therefore strongly decoupled from the small grain disk parameters. 

Having verified that we can reproduce the All result, that the All inferred gap radius and 
depletion factor for large sub-mm emitting grains is a good match to the models produced by 
both codes, and that this result is unaffected by variation in small grain parameters, we treat the 
All large grain disk parameters (with the exception of disk inclination, for which we adopt our 
measurement) given in Table 2 as fixed throughout the remainder of our modeling efforts. 
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4.3. SED Modeling 

There are sufficient data in the literature to construct a well-sampled observed SED for SR21. 
In addition to 2MASS J, H and K photometry, SR21 photometry has been obtained with Spitzer's 
Infrared Array Camera (IRAC) and Multiband Imaging Photometer (MIPS) at 3.6, 4.5, 5.8, 8 
and 70/xm, the AKARI spacecraft at 18/im, the Infrared Astronomical Satellite (IRAS) at 12, 25 
and 60^m, the Submillimeter Array (SMA) at 880//m, and the James Clerk Maxwell Telescope's 
Submillimetre Common-User Bolometer Array (SCUBA) at 350, 450, 850 and 1300/im. 10-35;um 
spectra from the Spitzer Infrared Spectrograph (IRS) and 52-97//m from the MIPS SED mode are 
also available. These are shown overplotted on all SEDs generated by our models. Although we 
show it for completeness, the MIPS 70/um photometry point has been excluded from our fits due 
to its inconsistency with both the MIPS SED spectrum and IRAS 60/im photometry. All reported 
photometric points and spectra have been extincted according to a standard ry=3.1 interstellar 



extinction law (Cardelli et al. 1989). The extinction for SR21 has been estimated to be as high as 
Ay=9.0 ( Prato et al.||2003 ), but All have argued that the published photometry is most consistent 
with a lower value of A\/~6.3. We also find that lower extinctions are most consistent with the 
observed optical/NIR SED, so we adopt the same Av=6.3 value for SR21. 

The SR21 SED shows circumstellar dust excess emission longward of the stellar blackbody 
peak, as evidenced by the long wavelength excess (see Figure 12). Unlike protoplanetary systems, 
which tend to have flatter transitions between the stellar and dust peaks in the SED indicating 
significant dust populations at all temperatures, the SR21 SED has a distinct double peaked profile, 
reaching a local minimum near 12//m. This suggests that there is a dearth of hot grains in the disk, 
a conclusion supported by the observed sub-mm cavity. Debris disk systems, on the other hand, 
are much more depleted of small grains, and their emission traces the stellar photosphere out to 
wavelengths of 20/um or more. In contrast, the SR21 SED rises rapidly beyond 12;um, indicating a 
healthy population of warm grains that emit at A>12/im. 

Shortward of the 12/im local minimum in the SR21 SED, the emission is also greater than 
photospheric, which indicates a small population of very warm grains (a hot inner disk component) 
or, alternatively, a warm companion. The SR21 SED shows little evidence of a 10/im silicate 
feature, while many similar transitional disk systems exhibit strong, broad silicate emission at this 
wavelength. A good disk model for SR21 therefore needs to reproduce all three of these distinct 
SED features: (a) the steep 12-20/zm rise, (b) the sub-12//m excess and (c) the lack of a strong, 
broad 10//m silicate feature. 

The Whitney code produces stable SEDs with relatively few (~10^) photons, making it much 
less computationally expensive than generating images/radial profiles. We were therefore able to 
explore much of the available parameter space for SR21 to find several classes of reasonable models 
capable of reproducing the observed features of the SED. 
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4-3.1. NIR Excess Fits 

The NIR excess shortward of 10//m is a distinct feature in the observed SED. The (small) excess 
in this region was originally attributed to a 0.07-0.25AU warm inner disk component. Indeed, 
adding this component to our models produces a well-fitting NIR excess. However, the MIR 



observations of Pontoppidan et al. (2008) are somewhat inconsistent with the presence of this 
interior hot disk component. Their data show a sharp truncation in gaseous CO at 7AU, suggesting 
a high flux of UV photons. This feature would be difficult to produce with an inner hot disk 
component at r<7AU capable of intercepting them. Furthermore, allowing the large grains to fill 
the cavity in to 0.07AU in our models, even while depleting them by a factor of 10~^, results in 
point-like sub-mm emission at the disk center, which is not present in the All observations. 

The presence of the warm companion predicted by E09 is another potential mechanism for 
creating excess NIR flux in SR21, and in this case allows for a larger inner small grain disk truncation 
radius. Our observations can neither confirm nor exclude the presence of a companion, yet they do 
provide indirect support of its existence, as detailed in Section 3.3. 

Otherwise identical models with an (a) 7AU truncation and a warm (T=730K, R=30rQ, E09) 
blackbody companion (b) a disk that extends inward to the sublimation radius at 0.07AU are shown 
in Figure 12 to demonstrate that both assumptions reproduce the observed NIR excess. Without 
any means of probing disk structures interior to the saturation radius of O'.'l in our observations, 
we cannot break the degeneracy between inner disk extent and the existence of a companion in the 
NIR SED fits using scattered light imagery. 

We hold a slight preference for the warm companion hypothesis based on the P08 CO results 
and our own observations of excess intensity and deviation from polarization axisymmetry (indi- 
rectly supporting the companion hypothesis and described in detail in Section 3.3), so we assume 
a 7AU truncation in all of our scattered light models. We add the warm companion component to 
all of our modeled SEDs in order to reproduce the observed NIR excess, but acknowledge that a 
significant inner disk extent is also a plausible mechanism for creating this feature. We do not add 
polarized intensity emission from the companion into our simulated images because the polarized 
intensity contribution of such a companion is (a) not well constrained by models and (b) should 
not have a significant effect on the azimuthally- averaged radial profile. 



4-3.2. lO/im Silicate Fits 



The final feature of note in the SR21 SED is the lack of a prominent 10/im amorphous silicate 
feature, which tends to be rather strong and broad in Herbig disks and is believed to originate from 
the surface layers of the inner (i20AU) disk ( |van Boekel et al.||2005"| |Meeus et ar]|2001| ). The IRS 
spectrum spans this region of the SED and also that of the weaker 18/im silicate feature without 
any evidence of a strong line flux at either wavelength. 
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We find that the majority of the 10/im siUcate emission comes from the inner disk wall in our 
models, and we can suppress it by either restricting the scale height of the small grain inner disk 
wall or by restricting small grains to a very small proportion of the total mass (<1%). Generally 
speaking, many otherwise well-fitting models overproduce 10//m silicate emission. Changing the 
small grain prescription may be warranted in future modeling attempts. Specifically, the silicate 
features could be suppressed by either increasing the grain sizes or limiting the proportion of silicate 
grains in the disk. We leave this for future modeling studies. 

Our scattered light imagery cannot give any information about the size of the small grain 
disk wall where the lO/xm feature originates, whether it lies at 7AU or 0.07AU (the sublimation 
radius). We limit the size of the inner disk wall to small scale heights in all of our scattered light 
models in order to suppress the 10/im feature, but in some cases not severely enough, resulting 
in overproduction of 10//m emission in the SEDs. We could presumably correct this by running 
additional models with smaller inner wall scale heights or reduced small grain populations, but 
have elected not to focus on reproducing this feature precisely for models that are poor fits to the 
radial polarized intensity profile due to the computational expense. The quality of fit to the lO^m 
feature was taken into account in choosing the final best-fitting disk model, described in section 
4.4 and shown in figures 15-16. 



4.3.3. Fits to 10-20fim rise 



The steepness of the rise from 12-25/im is perhaps the most unique and significant feature of 
the SR21 SED and is the primary feature that led to its "transitional disk" classification. The 
presence of this feature requires that there be a sharp increase in disk emitting area (e.g. a 
"wall") at intermediate (20-40AU) radii. However, the precise geometry of the "wall", as well as 
its composition, is not well constrained by the observed SED. 

The simplest way to create a "wall" is to allow for depletion inside of the r=36AU gap, so 
that the disk makes a transition from optically thin to optically thick at this radius. We know that 
the the disk is optically thin at r<36AU for long wavelengths (and therefore in large, thermally- 
emitting grains), and optically thick beyond; therefore, the large grains certainly do present a "wall" 
at 36AU. All of the models that we consider are therefore "gapped" in large grains. 

However, the sub-mm data places no such constraint on small grain depletion. In fact, there 



is reason to believe that small grain dust filtration occurs at gap radii ( Zhu et al. 2012 ; Rice et al. 



2006), and that material must survive into observed sub-mm dust gaps, as a number of gapped 



disks are still observed to be actively accreting ( Najita et al.|2007 Espaillat et al.||2012 Salyk et al. 



2009). The observed H-band radial profile of the SR21 disk would seem to support this hypothesis. 



as it shows no evidence of a change in small grain disk properties at r=36AU. Until our Subaru H- 
band data were obtained, there were no spatially resolved data at shorter wavelengths to constrain 
the properties of the small grains and the amount to which they are depleted inside of the sub-mm 
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cavity. 

The simplest hypothesis, put forward by previous SED modelers (Brown et al. 2007, All, 
E09) is that the small grain disk follows the same structure as the large grain disk, including the 
r<36AU gap. This hypothesis reproduces the observed SED, but it is inconsistent with our new 
Subaru H-band polarized intensity profile. Simulated radial profiles and SEDs for various small 
grain depletion factors are shown in Figure 13. They show that although small grain depletion 
in the r<36AU region reproduces the 12-20^m rise in the SED, it is strongly inconsistent with 
the observed H-band radial profile. All models in which the inner disk is strongly depleted, and 
therefore optically thin, suggest that we should have seen an excess in H-band emission at the 
cavity wall and a turnover in the radial profile at this radius. 

Undepleted and low depletion models, on the other hand, are inconsistent with the observed 
SED. Additionally, their radial profiles are too shallow to be consistent with the H-band data. 
Appendix A details our efforts to find a scenario in which the small grain disk, whether optically 
thin or optically thick at r<36AU, contributes to the 12-20/im rise in the SED. 

In the end we find that the steep 12-20/im rise can be created with a disk depleted only in 
large grains in its interior. Our data suggest that the small grain disk is not "gapped" . 



4.4. Scattered Light Imagery 

At first glance, the H-band radial polarized intensity profile is difficult to reconcile with the 
large, heavily depleted sub-mm cavity, as it spans the region of the sub-mm wall without any 
evident break or turnover. In fact, the principal features of the observed radial polarized intensity 
profile, namely its smoothness and its steepness (~r~^), were not reproducible under any of the 
scenarios we investigated in which the small grain disk contained a depleted cavity. Furthermore, 
we were unable to reproduce the radial profile under any scenario in which the H-band scattered 
light originates from the surface of an optically thick small grain component, whether gapped or 
not. This effort is described in detail in Appendix A, and representative models for each variety of 
optically thick disk model investigated are shown in Figures 17-19. This is a surprising result, as 
we had no reason to expect that transitional disks are optically thin at short wavelengths. 

We find that the steepness of the radial profile requires that we invoke an optically thin small 
grain component that spans the entire disk. In order to create this, we must either (a) relax 
the All assumption that the small grains compose 15% of the disk's dust mass and restrict their 
contribution to <1%, making the small grain dust disk optically thin throughout, or (b) invoke a 
third optically thin disk atmosphere or envelope component that has a larger vertical extent than 
either the large or small grain disks. 

Under scenario (a) the optically thick "wall" producing the 12-20/im SED rise must be com- 
posed entirely of large grains in order to allow the small grain disk to be optically thin. As we 
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noted in section 4.3, we do not require a small grain wall in order to reproduce the 12-20/xm rise 
in the SED. In fact, we can reproduce the SED structure at A>12^m with a large grain gap alone. 
We do need to increase the scale height of the gap wall in order to provide enough flux at these 
wavelengths (relative to models where an optically thick small grain component contributes to the 
rise), but not by such a degree as to lead us to declare these models unphysical {Zwaii ^lAU). 

Since an optically thin small grain disk means that we inevitably "see" the large grains at 
the midplane, these models require that the large grain disk not scatter efficiently at H-band. 
Otherwise, they suffer from the same local maximum problem that led us to exclude "Gapped" 
small grain disks. As can be seen in Figure 14, which shows a scenario under which the small 



grain disk is optically thin throughout, the Wood et al. (2002) dust prescription chosen for the 



large grain disk (and found to be a good match to that of All) does scatter weakly at H-band, 
resulting again in a turnover in the radial profile at 36AU. These models would require an alternate 
grain size distribution and/or composition in order to prevent a local excess at 36AU. We have not 
endeavored to explore this possiblity here. 

We disfavor explanation (a) because it requires that small grains be optically thin throughout 
the vertical extent of the disk. An entirely optically think small grain disk disk is very unlikely 
given the small amount of disk mass allowed and the observed optically thick r>36AU large grain 
component. In order to be optically thin throughout, we require not only that <0.1% of the disk 
mass lie in the small grain disk, representing a likely unphysical degree of grain growth or small 
grain depletion, but also that the large grain disk's grain size distribution contain no grains that 
scatter at H-band, an equally unlikely scenario. 

Under scenario (b), which we believe to be far more likely, an envelope component lies above 
the optically thick small grain dust disk with a significant vertical extent. This envelope contains 



very little total mass (on the order of 10 Mq) and is composed of small ISM-like grains (Kim 



et al.|1994 ) Representative model outputs for this scenario are shown in Figure 15, and a schematic 



of the disk structure is shown in Figure 16. 

Although we call this third disk component an "envelope" , it is not clear that it fits into the 
canonical definition. There has been no evidence in the literature to suggest that SR21 is still 



actively accreting, and the only constraints on the accretion rate are upper limits (Natta et al 



2006). This is in contrast to a number of other gapped disks that do show evidence of accretion 



dSalyk et al.|[2009l ). 



It seems that this optically thin component may be more readily explained as a vertically 
extended disk "atmosphere" of sorts, rather than an accretion envelope. Whether this is part of 
the natural structure of the disk, a remnant of the envelope, or a disk wind launched from the 
optically thick disk surface cannot be determined from our data. In fact, our observations are 
relatively insensitive to its precise geometry. The radial profile of the optically thin component 
follows an power law in all cases where the disk a-(3=-l, which allows for a large range in 
flaring and density structures. 
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Under this three component disk model, the radial profile is still sensitive to the structure of 
the underlying optically thick small grain disk surface. In particular, the addition of the optically 
thin component that gives rise to the power law cannot mask a discontinuous small grain 
distribution at the midplane, even when the optically thin component is allowed to contribute an 
order of magnitude more scattered light than the optically thick component. A large number of 
optically thick disk models were investigated in this effort, and are described in detail in Appendix 
A. 



5. Modeling Discussion 

Our requirements to fit the radial profile are threefold. First, wc require that there be enough 
small grain material inside the 36AU large grain cavity to scatter direct NIR starlight before it 
reaches the cavity wall, effectively shadowing it. Although this gap must be optically thin to sub- 
mm light in order to reproduce the SMA observations of All, we find no viable scenarios in which 
it can be optically thin at shorter wavelengths and remain consistent with the observed H-band 
radial polarized intensity profile. 

Secondly, the smoothness of the observed H-band polarized intensity profile suggests that the 
surface of this underlying optically thick small grain disk must follow a smooth surface geometry. We 
cannot fit our data under scenarios where there is a marked effect, either in differential depletion or a 
scale height discontinuity, of the large grain dust gap on the small grain distribution. This is perhaps 
the most interesting and puzzling result of our observations and analysis - that the small H-band 
scattering grains in the SR21 disk must follow a smooth distribution that is relatively unaffected by 
the goings-on at the disk midplane. This supports scenarios under which dust filtration occurs at 
the large grain gap wall, allowing the small grains to be carried into the gap along with the disk's 
gas. 

Finally, we find that the steepness of the H-band radial profile of SR21 is indicative of an 
optically thin small grain disk component. In order to remain consistent with the first two criteria, 
this requires that we invoke a third vertically extended optically thin disk component overlying the 
optically thick portion of the disk. 

As noted in Section 2, the AO system performance was sub-optimal during the SR21 obser- 
vations, and uncorrected effects such as those of a polarized seeing halo may mimic disk emission, 
particularly close to the star. Thus it is possible that the true SR21 radial profile flattens out in 
the interior and we are mistaking uncorrected PSF effects for disk emission in the inner region. 
However, all pixels in the saturated region were excluded from the radial profiles, and the PSF 
artifacts do not extend beyond the innermost bin in our radial profile plots, so this is unlikely. 

The computational expense of scattered light models means that we did not investigate the 
scattered light parameter space as systematically as the SED parameter space, and a more thorough 
and systematic exploration of the available parameter space is warranted to ensure that our solution 
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is unique. 

It could be argued, however, that ADI observations with more field rotation, higher AO per- 
formance (Strehl Ratio) H-band observations and/or further multiwavelength data that resolve the 
mysterious SR21 "gap" for other varieties of disk grain are equally warranted. Observing methods 
that are sensitive to sharp disk features (e.g., ADI) could help verify whether there is any local 
excess in scattered light at the 36AU dust wall in addition to resolving the potential companion. 

It is also worth noting that all of our disk models have assumed axisymmetry. Unresolved 
spiral structure, emission from a companion, disk warps, etc. could be present in the SR21 disk 
and may have an effect on the radial profile that has not been examined in this analysis. 

Uniformity of grain composition as a function of radius is another aspect of the models that is 
perhaps unphysical. Compositional differences, such as higher volatile abundances in the gapped 
region, could allow for a more marked effect of density on surface brightness in the inner disk, 
allowing optically thick large-/? large-a models to increase more steeply towards the star. 

6. Conclusion 

We have spatially resolved the SR21 transitional disk for the first time at NIR wavelengths. 
Our data show the signature of an optically thin small grain disk component, and place strong 
constraints on the underlying distribution of any optically thick surface grains. 

Our Subaru data resolve the disk in scattered polarized light for stellocentric r^C/.'! (>12 AU). 
Extended polarized emission is present above noise level out to r^^'.Q (~ 80 AU) and exhibits a 
centrosymmetric pattern of polarization vectors. The azimuthally-averaged radial profile of the 
polarized intensity emission reveals that the disk flux drops off steeply (r~^) and smoothly with 
radius. The disk is also evident in the PSF subtracted total intensity image, and the range of 
viable scale factors for PSF subtraction give us an estimated polarization percentage for the disk 
of 10%<P<27%. 

We confirm with an independent model that scenarios in which large midplane grains are 
depleted by a factor of ~10^^ inside of a '-^36 AU gap reproduce the SMA observations of All; 
however, small grains cannot be depleted to the same degree inside the cavity. 

Our analysis eliminates a variety of scenarios under which there is a discontinuity in the small 
grain disk at the 36 AU large grain gap radius, both in terms of physical likelihood and ability to fit 
the observational data. In order to be consistent with the lack of an equivalent gap in the Subaru 
HiCIAO H-band polarized intensity images and radial intensity profiles, there must be a decoupling 
of large and small grain disk parameters in the SR21 disk. Our data rule out scenarios in which 
the SR21 disk contains an optically thin small grain dust gap with a wall at 36AU. 

Although we can neither confirm nor deny its existence definitively, we present several pieces 
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of anecdotal evidence in favor of a companion in the SR21 disk. These include an intensity excess 
near the appropriate location of E09's hypothesized companion, a deviation from centrosymmetry 
in the polarization vectors in this region, and an apparent migration of the disk major axis from 
inner to outer disk, suggesting a disk warp or potentially unresolved spiral structure. 

We conclude that the best fitting models for the H-band radial polarized intensity profile, the 
observed sub- mm emission and the multiwavelength SED for the SR21 disk invoke three separate 
components. These axisymmetric, uniform radial grain composition scenarios suggest an optically 
thin envelope or atmosphere component lying above an optically thick, geometrically smooth small 
grain disk component. Further constraints on the geometry of this very interesting transitional disk 
may be obtained by follow-up at higher spatial resolutions and at other wavelengths. 
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Table 1: Position Angle Esti mates 



13 A 


Radius probed 


Reference 





ATT 




15±4 


7.0±0.4 


P08 (v=l-0) 


16±3 


7.6 ±0.4 


P08 (v=2-l) 


44±7 


llzbl 


E09 


61±8 


15±1 


E09 


76±8 


12-25 


this work 


96±7 


25-65 


this work 


100 


30-120 


All 


105 


30-120 


Brown 2009 
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Model Input Parameters 



Parameter 

^star 

Ay 

fgap 

i 

rc 



Model Input (s) 
Input Assumptions 
2M0* 
3.15R0* 
5750K+ 

6.3 
O.OO6M0* 
none, 18, 36* AU 
rs„6-7AU 
200AU 
14, 22* 
15AU* 



composition 
Zioo 

Oiouter 
1^ outer 

Oigap 

Pgap 
5 

composition 
id 
Zioo 

O^outer 
Pouter 

a 



gap 
figap 
6 



Large Grain Disk Parameters 



Wood et al, 



2002 ^- 



0.85* 
1.52AU* 
2.15+ 
1.15+ 
2.15+ 
1.15+ 

lo-^lo-^lo-'^^,lo-^ 10"^, 10-3, 10-2, lo-^ i 

Small Grain Disk Parameters 



Kim et al. 


1994 


Cotera et al. 


2001 



io-^io-sio-6,io- 



0.15 
2-15AU 
0.5-4 
1-3.15 
-5-10 
1-3.15 

10"^, 10-3, 10-2, 10-1, 1, 2, 5, 10 



Table 2: The range of Whitney model input parameters investigated for SR21. Where several 
discreet values were attempted, the table entry is given in the form of a list, and for parameters 
where a range of values was attempted, the table entry gives the full range of attempted values. 
Parameters identical to those of All are marked with an asterisk, while equivalent parameters 
(rounded for convenience, compatibility, or translated to Whitney code syntax) are marked with a 
+ sign. Rc is the characteristic radius for the exponential density cutoff as defined in All, f^ is 
the fraction of the disk's dust mass, Zioo is the disk scale height at 100 AU, a is the exponent for 
the radial midplane density dependence, /3 is the exponent for the radial scale height dependence, 
and 5 is the depletion factor of the gapped region relative to the outer disk 
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SR21 Polarized Intensity 




Fig. 1. — Median polarized intensity (PI) image of SR21 at H-band. The gap radius inferred by 
All and the companion orbit inferred by E09 are indicated with dashed elhpses. These are circles 
inclined at 14° and rotated to a major axis PA of 100° (All) and 60° (E09) (indicated with the 
solid white and red lines). North is up; east is left. The central O'.'O? were saturated in the raw 
images and have been masked out. A diffraction limited spot is shown at the bottom right for 
reference. Modeling suggests that the emission deficit ~0'.'2 south of the central source is a result 
of the disk geometry and is not a physical structure. 
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Fig. 2. — Left: Polarization vectors overplotted on the median polarized intensity (PI) map of 
SR21. A general centrosymmetric geometry is evident in the SR21 vectors, suggesting that this 
emission does indeed trace scattered light from a circumstellar disk. Right: Polarization vectors for 
the PSF star HD 148212 overplotted on its median PI image, which has been scaled by a factor of 3 
to compensate for the difference in exposure time. In this case, the vectors do not show a coherent 
pattern, and PI emission is weak and compact. This is interpreted as photon noise residuals and 
is used to estimate the error in the polarized intensity profiles shown in figure 3. In both cases, 
vectors are scaled by the strength of the polarized intensity signal in the region and are not plotted 
when they fall at or below noise level. Vectors whose lengths correspond to 500, 1000, 1500 and 
2000m Jy/arcsec^ are shown in the bottom right of the SR21 panel as a scale reference. 
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Fig. 3. — SR21 polarized intensity image (background) and polarization vectors (overplotted) after 
first order correction for a seeing limited halo component. This method reveals the polarization 
vectors to be more centrosymmetric, particularly along the disk minor axis) than originally inferred. 
Local deviations from centrosymmetry persist in the northeast and southwest quadrants of the disk, 
where unresolved disk structures and/or companions could induce local polarization structure that 
deviates from the larger centrosymmetric pattern. The process of subtracting this halo component 
has introduced additional PSF artifacts into the background PI map, most evident in the negative 
astigmatic pattern centered on the star. We have thus chosen to model the uncorrected polarized 
intensity map, which we believe more closely resembles the true structure of the SR21 disk. Vectors 
whose lengths correspond to 500, 1000, 1500 and 2000m Jy/arcsec^ are shown in the bottom right 
of the SR21 panel as a scale reference. 
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Fig. 4.— Intensity images of SR21(left) and HD 148212 (right). HD 148212 has been scaled 
according to differences in exposure time and intrinsic H-band magnitude (based on published 
2MASS values) to match SR21 for subtraction. Some uncorrected astigmatism is evident in the 
four lobes of the PSF in both stars. 



- 32 - 



SR21 Polarized Intensity Isophotes 
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Fig. 5. — Elliptical isophotal fits to the SR21 Polarized Intensity image. The fits are shown 
overplotted on the PI image in white starting at a semimajor axis of 10 pixels (well outside the 
saturation radius) to a semimajor axis of 50 pixels (where the PI fiux drops to noise level). The 
major axis of each isophote is overplotted as a straight white line. The average value of the major 
axis position angle is clocked by ~15° from the inner to outer disk, suggesting unresolved disk 
structure. The ellipticity of the isophotes was also used to estimate the inclination value of 14° 
(where 0° would be face-on). 



- 33 - 



1000 



SR 21 Radial Polarized Intensity Profiles 

|| I I I 

' SEEDS azimuthally averaged 

I SEEDS major axis 

I SEEDS PSF-subtracted 

I r^' profile 

I f^' profile 
36ALJ sub — mm cavity wal 

' sky background leve 



100 



A 



A' 



4\ 



A 



10 
0.1 



0.2 



0.5 



0.4 



0.5 



0.6 



radius in arcseconds 



Fig. 6. — Radial polarized intensity profiles for SR21 binned to the H-band diffraction limit (C.'S) 
and shown as triangles. The black filled triangles represent the azimuthally-averaged radial profile 
and the red open triangles represent the profile as measured in an 11 pixel wide region centered 
on the major axis of the disk. The blue open triangles represent the profile of the PSF-subtracted 
image. All three profiles follow a very similar slope, with the major axis profile systematically 
higher in flux per pixel by 20% relative to the azimuthally-averaged profile, and the PSF subtracted 
image systematically lower, as is to be expected. Because the images were saturated interior to 
(X'O?, radial profiles are only shown for bins that fall completely outside of the saturated region 
(from (X'l outward). Also shown is a canonical profile in green and a r^^ profile in blue. The 
profile is a better fit to the data, although there is some evidence that it becomes shallower 
in the outer disk beyond 0'.'45. The veracity of this shallower profile in the outer disk is brought 
into question by the steepness of the PSF subtracted profile, suggesting that the flattening in the 
outer profile may be due to halo effects. Errors were estimated from the photon noise in the four 
individual channels used to create the PI image summed in quadrature (black error bars). The 
location of the sub-mm cavity wall at 36AU and the noise floor are indicated with vertical and 
horizontal dashed lines, respectively. 
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SR2' PSF Subtracied Intensity 
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Fig. 7. — H-band PSF-subtracted intensity image with the 18AU orbit of the companion proposed 
by E09 overplotted. There is a distinct excess of emission in the Northeast quadrant (A) at 
a location that is consistent with their hypothesized companion, although its proximity to PSF 
artifacts makes robust identification impossible with this data. There is also a northward extension 
of the disk that is nominally present in the PI image (B). The image has been smoothed with a 
FWHM=6 pixel Gaussian in order to suppress artificial structures induced by speckle noise. 




Fig. 8. — Left: The 880;um SMA image of All. Right: A large grain disk model with All's best-fit 
input assumptions (5=10"^, M^j^^ = O.OO6M0, rgap=36AU, etc.) produced by the Whitney code. 
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6=10-6 
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Fig. 9. — Highly depleted ((5=10~^-10~^) large grain disk models produced by the Whitney Code. 
From left to right, the panels are (a) the observed SR21 SMA image of All (b) the final Whitney 
model image, where the raw model output was sampled with the same spatial frequencies observed 
with the SMA (c) the residuals and (d) the model fit (in red) to the SMA visibilities. Input 
parameters were identical to those of All, except that the depletion factor (6) for large grains 
inside of 36AU was varied. The depletion factor is given above each model row. A scale factor was 
also applied to each model image to equalize the model and SMA fluxes before subtraction from the 
observed image. This factor is given in panel (b) of each model. The minor discrepancies between 
modeled and observed fluxes can be easily explained by slight differences in input opacities. The 
bottom image is the same best fit model shown in Figure 7. 
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Fig. 10.— Moderately depleted ((5=10"5-10"3) large grain disk models produced by the Whitney 
Code. See Fig. 9 caption for details. 
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Fig. 11.— Minimally depleted ((5=10^2-10°) large grain disk models produced by the Whitney 
Code. See Fig. 9 caption for details. 
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Fig. 12. — Model NIR SEDs produced by otherwise identical models with and without (bottom and 
top panels respectively) the addition of a 700K blackbody companion to the SED and with either 
a 7AU inner disk truncation (right) or a disk extending inward to the sublimation radius (left). 
These reveal that the A <12^m NIR excess can be reproduced either by extending the disk inward 
to the sublimation radius (upper left) or by truncating the disk at a larger radius and adding in 
the NIR contribution of a TOOK blackbody companion (lower right). 



-40- 




Fig. 13. — Top: model SEDs of varying r<36AU depletion factors overplotted on the observed SR21 
SED. Five "Gapped" model runs of different depletion factors (5)are shown with solid colored lines, 
ranging from 10~^ depletion in red to zero depletion in purple. In the case of the heavily depleted 
6=W~^ and 5=10~^ models, the NIR excess in the SED can be reproduced either by extending 
the disk inward to the sublimation radius or by adding the contribution of a 700K blackbody 
companion. This is discussed in detail in section 4.3.1 Bottom: H-band radial polarized intensity 
profiles for the same models overplotted on the observed profile. The flux of the undepleted model 
was fixed to the data at r=36AU, and all other models scaled by the same amount, in order to 
highlight changes in the radial profile at the cavity wall. None of the models provide a good fit to 
the data, revealing that the small grain disk is not only undepleted, its geometry is different from 
that of the large grain disk. 
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Fig. 14. — Sample model outputs for a representative optically thin disk model. The lefthand panels 
represent the density of large grains (upper) and small grains (lower). The large grain model is 
identical to the All best-fit model described in Section 4.2 except that the dust wall at 36AU has 
been inflated by a factor of 2 in order to account for the lack of small grain contribution to the 
12-20/mMm rise in the SED. The small grain model has an inner truncation at 7AU. The righthand 
panels provide model fits to the observed SED (upper) and H-band radial polarized intensity profile 
(lower). In this case, the lack of 10//m silicate emission is due to the very small amount of small 
grain material populating the inner disk. As revealed by the radial profile, an optically thick small 
grain disk alone cannot mask the discontinuity at the 36AU large grain dust wall. As justified in 
the text and shown in Figure 15, an underlying optically thick small grain layer is needed to mask 
this discontinuity, while this optically thin layer is required to reproduce the r~^ radial polarized 
intensity profile. 
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Fig. 15. — Sample model outputs for a representative three component disk model. The lefthand 
panels represent the density of large grains (upper) and small grains (lower) . The large grain model 
is identical to the All best-fit model described in Section 4.2. The small grain disk consists of two 
components: (1) an underlying optically thick small grain disk of the 'Steeply Curved' variety with 
a 7AU inner truncation and (2) an overlying vertically-extended optically thin component that 
contributes thirty times more scattered light than the optically thick component. The righthand 
panels provide model fits to the observed SED (upper) and H-band radial polarized intensity profile 
(lower). The presence of the optically thick small grain disk acts to mask the discontinuity at the 
large grain dust wall that was revealed in Figure 14. In this case, the lOfim silicate feature is small 
due to the small vertical extent of the optically thick small grain disk wall and the small amount 
of material populating the more vertically extended optically thin small grain component 
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Fig. 16. — A schematic view of the three component disk model that provides the best fit to om' 
data. The vertical scale has been inflated by a factor of two relative to the horizontal. The origins 
of sub-mm light, H-band scattered light and the steep 12-20//m SED rise are indicated. 
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A. Optically Thick Scattered Light Models 

In this appendix, we detail our attempts to derive a disk model in which the H-band scattered 
light emission originates from the surface of an optically thick small grain disk. Although we 
ultimately invoked an optically thin component in order to explain our data, this was not the most 
obvious solution, and more than 100 optically thick small grain scenarios were investigated before 
we arrived at the conclusion that the steepness of the observer H-band radial intensity profile is 
due to the effect of an optically thin small grain disk component. 

Due to the relative expense of full scattered light models, we have not explored the entirety 
of the parameter space for the optically thick small grain disk solution. We therefore describe our 
efforts in this appendix, in the hope that they will help to inform future H-band scattered light 
modeling efforts. 



A.l. Gapped and Differentially Depleted Models 

If the r<36AU gap is optically thin at all wavelengths, as has been assumed by most previous 
models, then sufficient flux is incident on the dust wall at 36 AU to create a sharp rise in the SED 
from 12-20/im. We call the class of models with a depleted cavity that is optically thin in both 
large and small grains "Gapped" models. 

Our initial working assumption, though perhaps naive, was that the observed SR21 cavity 
would also be seen in small grains. As discussed in section 4.3, this hypothesis provides an excellent 
fit to the multiwavelength SED. However, if the cavity size and depletion factor suggested by the 
sub-mm data is carried over to the small grains that scatter at H-band, the simulated H-band 
polarized intensity radial profiles produced by the Whitney code (Figure 17) suggest that we should 
have seen an observable discontinuity in the H-band scattered light PI profile at the r=36AU sub- 
mm wall. We see neither a slope nor a flux discontinuity at this radius, and can immediately exclude 
the simple scenario in which the large and small grains follow precisely the same distribution. 

The simplest scenario that allows us to maintain a similar distribution for large and small 
grains is one in which the two populations are both depleted inside of the cavity, but each to 
a different degree (a different 6 value). Physically, this would seem to make sense, as the large 
grains concentrated at the midplane should be subject to different physical processes (perhaps 
most notably clearing by forming protoplanets) than the grains at the surface of the disk where 
H-band scattered light originates. Second generation small grains are also expected to be generated 



by the planet formation process (e.g., Zuckerman 2001), and so even if small grain depletion also 



occurs at the midplane, it could in some sense be compensated for. 

We began by exploring a variety of scenarios in which the small grains follow the same basic 
distribution as the large grains, are undepleted in the outer disk, but are depleted to a different 
degree than the large grains inside of the r<36AU cavity. Radial profiles for models following the 
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same basic distribution as the large grain disk (a=2.15, (3=1.15), but with depletion values ranging 
from 10~^ to zero (lO") depletion, are shown overplotted on the observed H-band radial profile and 
on the observed SED in Figure 13. 

In moderate gap depletion cases (5=10~^-10~^), the flux drops sharply inside of the gap and 
falls all the way into the saturation radius, which is inconsistent with our observation of a continually 
inwardly increasing profile. In cases where the depletion is high enough to make the disk optically 
thin ((5=10~^-10~*), the radial profile well inside the gap does increase sharply inward, however 
there is invariably a brightness discontinuity at the gap radius. This is due to the small amount 
of "wall" that is tipped into the line of sight by the disk inclination, which creates an excess at 
this radius that is evident in the radial profile. In low gap depletion cases (5=10*^-10"^), the radial 
profiles are relatively smooth, but too flat. 

The "Gapped" models provide excellent fits to the SED and sub-mm imagery. However, the 
increased emitting area of the gap "wall" invariably produces a local excess at C/.'S. In order for 
such models to fit the observed radial profile, two things would need to occur. First, scattered 
light emission from the wall (whose existence is required by the SED fits) needs to be somehow 
suppressed to avoid a sharp discontinuity in scattered brightness. If the disk were perfectly face-on, 
it may be possible to hide the excess emitting area of the wall, but the range of inclination estimates 
from a variety of multiwavelength work (All, this work, P08) suggest that this is unlikely. 

Second, they require that the optically thin grains inside of the gap and optically thick grains 
outside of the gap share a nearly identical radial PI profile despite very different physical conditions. 
Steep radial profiles can be produced by both optically thin and optically thick disks, but the 
apparent precise equality between the two makes this scenario seem somewhat contrived. We 
therefore exclude the "Gapped Disk" models as incompatible with the scattered light imagery 
except under the most contrived scenarios. 



A. 2. Scale Height Discontinuity Models 

We also consider scenarios that do not require the small grains to be depleted (and therefore 
optically thin) inside of the large grain gap. A step- like discontinuity in the small grain scale height 
could contribute to the steep SED rise, without requiring that the inner disk be depleted in small 
grains relative to the outer disk and allowing for identical (though discontinuous) distributions 
from inner to outer disk. We call this class of models "Scale Height Discontinuity" models, and 
the outputs for a representative model are shown in Figure 18. A related class of models utilize a 
localized "pufF' in the large grain outer disk wall. Infiated disk rims are often invoked to explain 



variable emission in disks (e.g., Flaherty et al. 2012). Under such models, the high incident flux 



on a disk wall causes higher temperatures at that location and "puffs up" the rim. If this were to 
occur at the 36AU large grain disk wall, it is conceivable that it would also inflate the small grain 
scale height at this radius. 
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At the limit of a perfectly vertical step-like discontinuity and a face on disk, there is no break 
in the radial profile between inner and outer disk at H-band because the optical depth of these 
grains is large and the wall is obscured by the disk surface as viewed from face-on. However, as 
with the gapped disk models, it is difficult to mimic this situation in an inclined disk, even one 
with low inclination like SR21. These models thus solve the problem of requiring a different disk 
geometry inside and outside the gap however they also invariably produce a local excess at the gap 
radius that is not seen in the data, just as the "Gapped" models did. 

In the case of puffed disk rims, the discontinuity that causes the puff is localized to the "wall" , 
so the disk soon settles back into its normal scale height distribution. These "puffs" may have 
a small radial extent, but again it is difficult to construct a situation in which they wouldn't 
be visible as an excess in the radial profile. On the other hand, the puffed disk rim shadows a 
portion of the disk immediately exterior to it from incoming stellar flux, where the area of the 
shadowed region depends on the height of the "pufF' . This local excess ( "pufF' ) followed by a local 
minimum (shadowed region) could be smoothed out in the polarized intensity radial profile if the 
total radial extent is much smaller than the diffraction limit at H-band, mimicking our observed 
smooth radial proffie. This precise balance again seems somewhat contrived, so we consider the 
puffed rim explanation to be relatively unlikely, and we do not include any of these models in our 
final fits. 

Owing to their continuous density distribution, brightness discontinuities in the "Scale Height 
Discontinuity" models are less drastic than those created by the "Gapped" disk models, so they are 
more likely to remain unresolved in our observational data. Non-diffraction limited performance, 
the peculiar shape of the HiCIAO PSF and the relative insensitivity of polarized intensity emission 
to sharp disk features (due mainly to low overall flux) could all conceivably smooth out such 
discontinuities in the radial profile to a greater extent than convolution with a 0.06" Gaussian PSF 
would suggest. However, we were unable to construct a disk model in which the radial profile 
was smooth enough to be consistent with our data and error bars, so we exclude these models as 
relatively unlikely candidates. 

A.3. Steeply Curved Models 

If we relax the assumption that the small grain disk needs to have precisely the same parameters 
inside and outside the gap, we find that there are cases where a turnover in the radial scale height 
exponent ^ at the 36AU wall can reproduce the steep SED rise. This situation does not require 
small grain depletion, nor does it require a discontinuity at the large grain dust wall. In such 
cases, the SED rise coincides with the region just inside of 36 AU where the disk is most steeply 
curved, and in that sense mimics a "wall". From a physical standpoint, this transition from a large 
to a smaller radial scale height exponent at the large grain dust wall would need to be explained 
somehow by the dearth of large grains at the disk midplane in this region or the physical mechanism 
causing disk clearing, which we do not endeavor to do here. We call this class of models "Steeply 
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Curved" models, and outputs for a representative model are shown in Figure 19. 

Under the third well-fitting SED scenario, a steep disk curvature (high /?) in the inner disk 
mimics a "wall" at the transition to a shallower outer disk profile. Without a physical discontinuity 
in the small grain surface at 36 AU, these models are also free of discontinuity in their radial 
polarized intensity profiles. However, all "Steeply Curved" models also produce flat (r^) radial 
profiles, which our observations argue against. This is because the geometric emitting area per 
pixel increases as the disk becomes increasingly flared, so pixels represent a larger amount of disk 
surface as you move outward. This mitigates the effect of incident stellar radiation, and makes 
the flux in the disk extremely flat. 

We find that the surface brightness of the disk is relatively insensitive to the radial density 
exponent a in our high-/3 models, even for very strongly inwardly peaked density profiles (e.g., 
r~^^). These models still show a flattening of the brightness profile in steeply curved regions that 
our data rule out. We therefore exclude these high-/? undepleted scenarios as poor fits to the H-band 
data. 

We find that the steepness of the radial profile significantly limits the disk geometry, particu- 
larly in the optically thick case. We tried more than 30 iterations of smooth optically thick disks 
with a wide range of a, (3 and a-f3 (the radial dependence of the surface density) values, and were 
unable to create a steep profile. 

The only optically thick disks we found that reproduce an radial polarized intensity profile 
are geometrically flat. As the geometrically flat case requires that the scale height of the disk be 
smaller than the radius of the star, it is impossible to create one in which the large grain disk 
required to reproduce the sub-mm flux does not stick up above the small grain disk. It is difficult 
to imagine a situation that would reproduce this, so we exclude these models on physical grounds. 



A. 4. Prospects for Envelope Component 

Since SR21 may still be weakly accreting (the upper limit of 10~^-^^ on the accretion rate is 
relatively high compared to other disks), we also explored the possibility that an envelope compo- 
nent is contributing to the flux in the interior of the gap. To provide observational constraints on 



the amount of circumstellar material, SR21 was observed on 2011 March 29 with SIRPOL (Kandori 
et al.|[2006 ), a JHKs-simultaneous imaging polarimeter that is mounted on the Infrared Survey Fa- 
cility (IRSF) 1.4m telescope in South Africa. We conducted aperture polarimetry on this data and 
the results are summarized in Table 3. Our data reveal total polarization levels that are consistent 



with the range of total polarization observed for YSOs ( Pereyra et al.|2009 ). We have not corrected 



the SR21 aperture polarization values reported in Table 3 for interstellar polarization, which could 
be signiflcant due to the high foreground extinction (Ay =6. 3). Moreover, because the observations 
were conducted in defocus mode to avoid saturation, the A and B components of SR21 (separated 
by 6") overlap in our aperture polarimetry, preventing us from using SR21B to determine how much 
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of the measured polarization is due to foreground dust and how much is local to SR21A and/or B. 
There are no other nearby objects in the field suitable for use as field star interstellar polarization 
probes. 

If we assume that the level of interstellar polarization and contamination from SR21B is min- 
imal, we can use the measured polarization Position Angles (PAs) at each band to probe of the 
optical depth of the disk, as optically thick disks show polarization PAs oriented perpendicular 
to the PA of the disk major axis, while the polarization PA and major axis PAs of optically thin 



disks are oriented parallel to one another (Whitney & Hartmann 1992 Pereyra et al. 2009). The 



disk major axis PA (86ibll°, this work; 100°, All) is roughly perpendicular to the measured po- 
larization orientations at all three bands, suggesting that the disk may be optically thick in the 
NIR. 

However, there are several caveats to be made. First, the 42° counterclockwise migration in the 
disk PA derived from isophotal fits (and described in section 3.1) suggest that an average PA is not 
a particularly good measure of the overall disk orientation. Secondly, these results depend on the 
robustness of our assumption that the interstellar polarization and contamination from SR21B is 
small. Finally, and most importantly, we have managed to eliminate all of the optically thick small 
grain dust disk scenarios that we have investigated thus far. In fact, we find that optically thin 
disks provide a much better fit to the H-band radial profile, as described in Section 4.4, although 
we cannot rule out, and in fact require, an additional contribution from an underlying optically 
thick small grain disk component. 



Table 3: SIRPOL Aperture Polarimetry 



Band 


Polarization 


PA 




(%) 


O 


J 


2.32± 0.17 


21.5±2.1 


H 


l.OlibO.05 


18.0±1.3 


Ks 


0.52± 0.10 


179.2±5.3 
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Sample Gapped, Smooth Model 
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Fig. 17. — Sample model outputs for a representative gapped disk model. The lefthand panels 
represent the density of large grains (upper) and small grains (lower). The large grain model is 
identical to the All best-fit model described in Section 4.2. The small grain disk, like the large 
grain disk, is depleted by a factor of 10^^ interior to 36AU, and it truncated at 7AU. The righthand 
panels provide model fits to the observed SED (upper) and H-band radial polarized intensity profile 
(lower). Although the small grain surface layer is continuous under this model, the heavy small 
grain depletion interior to 36AU allows a significant amount of flux to be indecent on the wall, 
creating a discontinuity in the H-band polarized intensity profile. 
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Sample Scale Height Discontinuity Model 

Large grain density distribution Spectral Energy Distribution 
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Fig. 18. — Sample model outputs for a representative undepleted small grain disk with a scale height 
discontinuity at the 36AU large grain dust wall. The lefthand panels represent the density of large 
grains (upper) and small grains (lower). The large grain model is identical to the All best-fit 
model described in Section 4.2. The small grain disk undepleted, with a scale height discontinuity 
at 36AU and an inner truncation at 7AU. The righthand panels provide model fits to the observed 
SED (upper) and H-band radial polarized intensity profile (lower) . The large amount of small grain 
material in the inner disk overproduces the 10/im flux, but provides an otherwise good fit to the 
SED. In this case the discontinuity in the radial profile at 36 AU is smaller than in the gapped 
models due to the relatively smaller amount of direct stellar flux incident on the small grain dust 
wall. However, the steep r~^ profile is not reproduced in this case, nor was it with any of the other 
radial density (a) and scale height (/3) exponents investigated for such scale height discontinuity 
models. 
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Sample Steeply Curved Model 

Large grain density distribution Spectral Energy Distribution 
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Fig. 19. — Sample model outputs for a representative model in which the region of the sub-mm 
cavity is more steeply curved in the small grain disk. The lefthand panels represent the density of 
large grains (upper) and small grains (lower). The large grain model is identical to the All best-fit 
model described in Section 4.2. The small grain model is undepleted interior to 36 AU and truncated 
at 7AU. The scale height exponent /3 was increased from 1.15 in the outer disk to 2.5 interior to 
36AU, allowing for a steeper curvature. The righthand panels provide model fits to the observed 
SED (upper) and H-band radial polarized intensity profile (lower). This small grain distribution 
completely masks any discontinuity at the 36AU dust wall created by an H-band contribution from 
the large grain dust prescription, however the steep inner disk carries a flat radial profile, which is 
inconsistent with the observed r~^ profile. 
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